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ABSTRACT 

This paper represents the Initial quantitative comparison of Inlet suppres- 
or far-fleld directivity suppression with that predicted using an appro;<lmate 
liner design and evaluation method based upon mode cutoff ratio. The experimental 
data was obtained using a series of cylindrical point-rnctinp inlet liners on an 
Avco-’Lycomi ng YF102 engine. The theoretical prediction program Is based upon sim- 
plified sound propagation concepts derived from exact calculations. These indicate 
that all of the controlling phenomenon can be approximately correlated with mode 
cutoff ratio which itself is intimately related to the angles of propagation within 
the duct. The objective of the theory-data comparisons Is to point out possible de- 
ficiencies in the approximate theory which may be corrected. After all theoretical 
<ri refinements liave been made, then empirical corrections can be applied, 
n 

INTRODUCTION 

There is a need for an accurate but simplified method for designing and evalu- 
ating supnressors without splitters for aircraft engine inlets. Methods have been 
published’’^ which can be used to estimate acoustic power reduction in exhaust ducts 
or in inlets with splitter rings. However, what Is really needed is a comprehensive, 
generally available, design or evaluation program which provides acoustic suppres- 
sion information in the far-field. 

Comprehensive suppressor evaluation is a very complex problem involving source 
acoustic oower description, suppressor effectiveness, modal scattering at several 
places in the propagation path and, finally, radiation to the far-fieid. The com- 
plete coupled system may never be amenable to analysis, except perhaps by some as 
yet undeveloped numerical analysis program. The method used in this paper represents 
a first order, uncoupled, analytical approach to each element of the propagation and 
radiation problem. Each element was found to be dependent upon the mode cutoff ratio, 
which provides a great simplification to the entire problem. Thin dependence upon 
cutoff ratio has been reoorted for maximum possible attenuation^, optimum impedance^, 
duct termination reflection^, far-field rr.diation directivity^ and modal density 
functlon^^. 

The purpose of this paper is to compare calculated noise suppression results 
using the simplified analysis, in its present state of development, with experimental 
Inlet suppressor data reported in reference 7. The comparison will be made on the 
basis of far-field attenuation directivity. This comparison, made on an angle-by- 
angle basis, is extremely demanding since it encompasses all elements of the problem 
starting with the description of the noise source and then proceeding with the sup- 
pressor performance and, finally, the sound radiation from the duct inlet. The far- 
field attenuation directivity is, however, the essential information needed for sup- 
pressor performance. 

’’Head, Acoustics Section; Member AIAA 
■^A^^osp^ce Engineer; Member AIAA 
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The future work concerning this simplified liner evaluation approach 
will proceed In two directions. The theory-data comparisons will be 
used to point out where further analytical effort Is required to encom- 
pass more of the essential physics Into the program. Also, an effort 
will be made to account for missing phenomena by empirical corrections 
to the theory. 
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three Inlet liner constructions (Table i) 

blade passage frequency, Hz 

liner backing depth, m 

speed of sound, m/sec 

modal density function, see equation (7) 

circular duct diameter, m 

orifice diameter in perforated plate, m 

sound attenuation, decibels 

maximum possible sound power attenuation, decibels 
sound attenuation using liner B3 only, decibels 
sound attenuation using both liners B3 and B2, decibels 
sound attenuation using liners B3, B2 and B1 , decibels 
quantity In circular contour correlation, see equation (A5) 
frequency, Hz 

acoustic Intensity Incident upon duct termination, N/MS 

acoustic Intensity transmitted out from termination, N/MS 
« 

length of acoustic treatment, m 

axial steady flow Mach number, free-stream uniform value 


spinning mode lobe number (circumferential order) 
number of modes 

exponent in cutoff ratio biasing function, see eiquations 
(6) and (8) 

far-field acoustic pressure amplitude, N/M^ 

far-fieid acoustic pressure amplitude due to multiple modes, 
N/M^ 

eigenvalue for circular hardwalied duct or absolute value of 
complex eigenvalue for softwail duct 

radius in circular damping contour correlation, see equation 

m 

circular duct radius, m 
perforated plate thickness, m 
inlet 1 iner (Table I ) 
radial eigenvalue 

cutoff ratio biasing function, see equation (8) 

AdB^/AdB 

l/7th power boundary layer thickness, m 
displacement boundary layer thickness, m 

6/ro 

frequency parameter, q = fD/C 
specific acoustic resistance 

resistance coordinate for center of damping contour 
optimum specific acoustic resistance 

optimum specific acoustic resistance without boundary layer 
(slip flow) 


radial mode order Index 
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a 

X 

Xc 

♦p 

fo 


mode cutoff ratio 

liner open area ratio 

phase of complex eigenvalue, degrees 

specific acoustic reactance 

reactance cooi i.!iiiate for center of damp! ig contour 
optimum specific acoustic reactance 

far-fleld angle measured from Inlet duct axis, degrees 
i)j at peak of principal lobe of radiation 
directivity angle without flow, degrees 

SIMPLIFIED APPWXIMATE SUPPRESSOR THEORY 


In this section the elements of the simplified suppressor theory will 
be outlined. Since this paper represents an Initial theory-data compari- 
son and not a final recommended procedure, many of the correlations used 
will be only referenced. Only the recent changes in the procedure and ex- 
pressions which are useful in later discussions will be listed explicitly. 
The subjects covered include suppressor performance, far-fleld radiation 
and duct termination effects, source power distribution and assumptions 
implicit in using the abo^e as uncoupled elements in a system which is ac- 
tual ly coupled. 


Suppressor Performance 

The sound attenuation obtained in a cylindrical suppressor is obtained 
by using correlations of the optimum impedance and maximum possible atten- 
uation. These feed into a circular approximation correlation of the actual 
damping contours which are used to estimate off-optimum suppressor perfor- 
mance. 

The concepts of optimum impedance and maximum possible attenuation are 
illustrated In Figure I. Equal damping contours representing exact calcu- 
lations for a particular mode (m = 7, P = 0 are represented by the solid 
lines. As damping is increased, the contours shrink in size and the limit 
is represented by the optimum impedance. The damping at this limiting 
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Impedance Is called the maximum possible attenuation for this mode. 

Pushing to higher attenuations would result In a contour which can be 
related to the next higher radial mode (y « 2 In this case). In ref- 
erence 8 It was discovered that modes with similar cutoff ratios had 
similar optimum Impedances. This Idea was pursued In reference In 
which a wide range of parametric study of optimum Impedance was made. 

The mode cutoff ratio defined In the suppressor section as, 

( 1 ) 

R ^(1 -Mg) cos 2(}) 

was found to correlate the optimum Impedance quite accurately. The 
only modal Information required Is, thus, the cutoff ratio. The op- 
timum Impedance components are found using the procedure reported in 
reference k and are a function of the frequency parameter (q) , bound- 
ary layer thickness, steady flow Mach number and, of course, cutoff 
ratio. 


The maximum possible attenuation was also shown to be approximately 
correlated by cutoff ratio alone^. The expression used here is, 
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Equation (2) shows that maximum possible attenuation is a function only 
of cutoff ratio, Mach number and duct L/D. Al though optimum Impedance 
can be a strong function of boundary layer thi ckness^ >9 , it was shown 
in reference 9 that the maximum damping Is insensitive to boundary layer 
thickness. 


( 2 ) 


Off-optimum suppressor performance is evaluated by an approximate 
circular damping contour correlation similar to that developed in ref- 
erence 10. This approximate correlation Is compared to exact calcula- 
tions in Figure 1 for a well propagating mode (large ^) . The correla- 
tion of reference 10 was semi-empIrIcal and was based upon calculations 
made for well propagating modes. When theory-data comparisons of far- 
field radiation attenuation were made. It became obvious that the cir- 
cular contours were not sufficiently broad for cutoff ratios near unity. 
An evaluation of the exact equal damping contours near cutoff was made 
and the correlation was altered to reflect this near cutoff broadening. 
The new correlation is presented In the Appendix. An additional change 
from the correlation of reference 10 is also found in the Appendix. 

When optimum resistance is reduced due to boundary layer refraction, 
the entire family of contours is also compressed. The new correlation 
Is compared to exact attenuation contours for a near cutoff mode In 
Figure 2. Notice that for low attenuation (larger contours) the con- 
tours do not close as they did In Figure I for high cutoff ratio. 
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Par-Field Radiation 

The basis of the far-fleld radiation expression used for the calcula- 
tion reported here is^, 


p2 R/ 2 sin 'll \/l-l/g^ fsin[ffn(sin ili - l/g)l { ^ 
TT^n [l/C^ - sln2 


( 3 ) 


which has the principle lobe peak located approximately at the angle 
given by, 


sin i|i ^ l/£ 

P 

where the mode cutoff ratio (?) Is defined by, 


(M 




- 2J1 
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(5) 


Equation (3) represents an approximation to the exact equation for a 
flanged duct with equal acoustic power per mode''. The modification to 
the exact expression to yield equation (3) was made mainly to allow the 
expression to be only a function of cutoff ratio (^) by the suppression 
of modal information (m, a). The approximations Involve using the sinu- 

soidal approximation to the Bessel function and assuming that m/a « I. 
Equation (3) represents an adequate approximation to the exact flanged 
duct expression except for the lower radials of the high lobed (larger m) 
modes where some error occurs'^. 

The far-field radiation pattern for multimodal radiation can be ap- 
proximated by, 

» 2 cos ij; sln'^ i|; (6) 

which was derived in reference 5 by consideration of the amplitude and 
location of the principal lobe of equation (3) and using the modal 
dens I ty function® 

N dC ~ ^3 


(7) 


(8) 


7 

5 

and thf cutoff ratio biasing function , 
P 1/f," 


Equation (6) expresses the effect of the summation of ail of the princi- 
pal lobes of radiation for the propagating modes. Equation (7) allows 
the estimation of the relative number of modes within a desired cutoff 
ratio range, and equation (8) allows the biasing of the source distribu- 
tion towards cutoff away from equal acoustic power per mode. Note that 
when n “ 0 equal acoustic power per mode Is used and equation (6) re- 
duces to the simple expression 2 cos \p. Exact summations, using the 
flanged duct radiation equation, of all the propagating modes were made 
for n ® 0, 1, 2 and 3 In reference 12. Equation (6) was shown to be 
an extremely good approximation. 

In a static test case fc** an engine inlet, where there Is no flow 
in the far-fleld, It was found that the far-field radiation Is shifted 
toward the Inlet axis. The position of the peak of the principal lobe 
of radiation was found to be given by ^3^ 

-M + i/TiT? 

cos (9) 

1 - 


where now the cutoff ratio must be defined by, 




ir n 

R /Tm7 


(10) 


This radiation direction shift is caused by the change in axial propaga- 
tion angle for a duct mode due to the convective velocity effect '3, /\n 
example of this effect Is that at Mq = -OA, a mode near cutoff (?«>' 1) 
propagates at 66° from the inlet axis. At Mp = 0, a near cutoff mode 
has Its peak propagation at 90°. This Idea was further generalized In 
reference to Include the entire radiation pattern. The corrected 
radiation directivity including the convective effect can be expressed as, 

-H„ + l/l - 


I - 


COS 


( 11 ) 
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where i{(p Is the directivity angle without flow. This directivity shift 
toward the Inlet axis was found to be essential to explain the experi- 
mental attenuation data. Unfortunately, this shift was made only In the 
final attenuation calculations and has not yet been Integrated Into the 
computer program. For example, equation (6) must be altered by the con- 
vective shift to accommodate the principle lobes of radiation, and the 
shift must be applied also to the side lobes of equation (3) (see ref- 
erence 5 l^or current handling of side lobes). It Is expected that the 
Integration of the convective directivity shift will not greatly alter 
the results presented In this paper; however, this has not been checked. 

One problem occurs when the convective radiation shift is applied. 

Since the flanged duct radiation theory applies only to 90° for a no-flow 
situation, after the shift toward the inlet axis Is made the theory will 
only extend to about 66° (see example above). Thus, an unflanged radia- 
tion directivity expression Is needed for this program. These results 
are available In the 1 1 terature'5> 16, 17^ however, they are too unwieldy 
for routine use In a multimodal radiation calculation. More significantly, 
these equations have not yet been reduced to a form containing modal In- 
formation only In the form of cutoff ratio so that they can be Integrated 
Into the current program. Effort Is underway to obtain a simple, adequate 
approximation for this purpose. 


A duct termination loss Is also used In this program. This termina- 
tion effect was derived In reference 5 and Is given fay, 


*T ^ k l/l- 1/g^ 
(1 + \/l-l/C^ 


( 12 ) 


Equation (12) contains only the mode cutoff ratio. It Is known to the 
author that a better expression containing also the frequency parameter 
(n) should be used, but Is has not been completed at this time. 


Source Power Distribution 

It Is known that the sound power attenuation provided by an acoustic 
liner Is strongly dependent upon the acoustic mode present In the llner^»^°. 
Likewise, with a multitude of modes present,, the attenuation should be a 
function of the distribution of acoustic power among the modes. The method 
presented In this paper does not use modes directly, but Instead bases all 
of the elements of the program upon mode cutoff ratio. Thus, modal power 
distribution Is replaced by cutoff ratio power distribution. This Is an 
Important simplification since although many modes may be present within 
a small Increment of cutoff ratio, It is not necessary to know the power 
of the individual modes but only the total power of all of these modes. 


’I 
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This simplification also has practical significance. Although experl- ’ 

mental measurement of modal amplitude and phase has been attempted for 

quite simple modal structures in rectangular'^ and annular^O ducts, It 

appears that If a large number of modes Is possible then direct modal 

measurement will be extremely difficult^®, if not Impossible. ! 

The only available method to estimate an acoustic power distribution | 

Is to observe the far-fleld directivity pattern and to Infer the source 
distribution which could have caused this radiation pattern as In ref- 
erences 5 and ll. This Indirect method Is usable for e single mode or 
for multimodal propagation where the modes are uncorrelated. It would 
be difficult to assess the radiation pattern of a few correlated modes 
but then, of course, direct modal measurement in the duct would be feas- 
ible. Equation (6) can be used for a comparison with the far-fleld ex- 
perimental directivity data and the exponent n can then be estimated. 

Equation (8) will then show the biasing away from equal energy per mode. 

This was done In reference 7 for the fan blade passage tone considered 
In this paper with the result that n 0, Another method that could 

be used when the principal lobes of radiation dominate the entire dlrec- 1 

tivlty pattern Is to compare directivity data at each angle to equation I 

(6). Since the far-field directivity would be dominated by radiation I 

from modes of a particular cutoff ratio range, then several discrete I 

corrections (away from equal power per mode) at several cutoff ratios ! 

could be obtained. This refinement was not attempted here due to the 
previously mentioned corrections which have yet to be developed for 
equation (6) to account for duct flow convective directivity shift. 

^ Thus, as a first approximacion, equal acoustic power per mode (n ■ 0) 
will be used for all of the calculations presented in this paper. 

- 

Inherent Assumptions of Present Program 

Several intrinsic assumptions are embodied In the methods as pres- 
ently reported. First, it Is assumed when many modes are present and 
that they are uncorrelated so the squares of the pressures (as In eq. 

(3)) can be summed (the squares of the sum of pressures must be used 
for correlated modes). This Is probably valid for randomly generated 
noise such as broadband noise or even tones generated by Inflow turbu- 
lence. 

Second, it is assumed that for a multitude of modes at all possible 
cutoff ratios, modal scattering has no net effect. For a single mode, 

the scattering occurring between a hard-walled and treated section could ! 

certainly affect the calculated suppression. Different radial boundary ] 

conditions are present between the hard- and soft-walled sections. Thus, 
differing axial and radial propagation angles will occur between the two 


i 



)0 


sections (see reference 21 for discussion of propagation angles In lined 
ducts). Therefore, although there Is no circumferential scattering, due 
to circumferential symmetry, there will, in general, be radial mode scat- 
tering occurring. This radial scattering Is probably a function of cut- 
off ratio, frequency parameter, Mach number and wall Impedance, Although 
scattering is known to occur, to simplify the problem It Is assumed that 
the scattering will balance out. For example, a -^Wen acoustic power 
versus cutoff ratio distribution may exist In the hard duct. For each 
group of modes near a given cutoff ratio there will be scattering to both 
higher and lower cutoff ratios. It Is assumed that as much power Is re- 
ceived from other cutoff ratio modes as Is lost to them. Thus, the power 
distribution in the soft duct (based upon cutoff ratio) Is the same as 
that In the preceding hard duct. A study of this phenomena must yet be 
made, and It Is possible that some of the discrepancies between theory 
and data, which are discussed In a later section, are due to neglect of 
this phenomenon. 

Reflections at impedance and area discontinuities are Ignored except 
for the crude duct termination effect given by equation (12). This equa- 
tion Is an approximation fo? the reflection of a given mode with the same 
mode Incident. Actually, there will be some radial mode scattering and 
reflection (similar to that discussed above in connection with the hard- 
soft wall Interface).. Subtle radiation pattern changes may occur which 
would lead to disagreement between theory and data. However, this effect 
Is Judged to be of ’ess importance than that at the hard-soft Interface. 

It is also assumed that the inlet Up shape has no effect upon the 
far-fleld radiation pattern. This effect Is probably small for Inlets 
with fairly abrupt lips and with high Mach number. The high Mach number 
causes the modal wave fronts to propagate at small to moderate angles 
to the axis’’ (onlv up to 66° for Mp « -O.k) and, put very simply, the 
emerging wave does not see an abrupt change in the lip. A gradual area 
change Is quite another case; the decreasing Mach number and Increasing 
radius can affect the propagation angles and, thus, the radiation angle. 
These effects are currently being studied and will be Incorporated Into 
the program when available. 


PROCEDURE FOR CALCULATING THEORETICAL ATTENUATION 


The procedure for combining the theoretical elements discussed In 
the previous sections will now be outlined. Figure 3 will be used to 
help Illustrate the calculation of the far-fleld directivity attenua- 
tion. The calculations are made for a particular frequency, duct Mach 
number, boundary layer thickness and duct geometry. The optimum im- 
pedance locus (the solid curve) is first generated by either exact 
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modal calculat* "s or more conveniently using the correlation equations In 
terms of cutoff ratio from reference h. The locations of the centers of 
equal mode number Intervals are found’ using the modal density function, 
equation (7)* Five of these locations are shown for Illustration In Fig- 
ure 3 but many more are used In the calculations shown later In this paper. 
Note that each of these points has a cutoff ratio associated with It and a 
preferred radiation direction Is noted for each as calculated from equation 
(l»). Plots of the experimental hard-wall far-fleld directivity are compared 
to equation (6) to get an estimate of the acoustic power distribution, Note 
that unless a bias (eq. (8)) Is applied, all points on the locus of Figure 
3 will contain equal acoustic power since they contain an equal number of 
modes, 

Each of the modal groups Is used to form a hard-wall far-fleld radia- 
tion pattern using the flanged duct equations of reference (5). This rep- 
resents a simplified Integration process over all modes whose principle 
lobes of radiation can contribute to a particular far-fleld angle, The 
procedure is the same as that used to derive equation (6). SIdelobes are 
also considered. The sidelobes of modes at all cutoff ratios can possibly 
contribute to the radiation at any given angle In the far-fleld. An Inte- 
gration procedure is also used for the sidelobes as presented In reference 
5. The hard-wall far-fleld radiation pattern Is, thus, established at the 
five angles represented In Figure 3> 

For the lined duct configuration, a wall Impedance of Interest (or an 
actual wall Impedance for a particular liner) Is chosen as In Figure 3. 

The damping of each of the five modal groups must then be calculated. The 
maximum possible attenuation Is calculated using equation (2) and, of course, 
the optimum impedance Is already known (the point on the locus), Off- 
optlmum attenuation (the circles and the rectangular point in Figure 3 do 
not generally coincide) Is then calculated using the equations in the Ap- 
pendix (Illustrated In Figures 1 and 2) defining the approximate circular 
damping contours. The original acoustic power at each cutoff ratio Is, 
thus, reduced by a certain amount by the liner and this reduced acoustic 
power is radiated to the far-fleld as In the hard-wall case. The difference 
between the hard-wall and soft-wall far-fleld directivities thus represents 
the far-fleld directivity attenuation. Equation (11) is then used to shift 
the directivity attenuations from the original angles ip^ to corrected 
angles \p. This accounts for the convective effect present with flow In 
the Inlet duct and no flow in the far-fleld (I.e., the static test case). 

The above discussion using five mode groups and, thus, five far-fleld 
angles was used for illustration only. The actual program uses differ- 
ent cutoff ratios to establish the function of suppression versus cutoff 
ratio. Interpolation is then used as needed. Far-fleld radiation patterns 
are calculated every 5 degrees from 5 degrees to 85 degrees. The convec- 
tive correction (equation (11)) must be applied to each of these angles. 
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DATA-THEQRY COMPARISONS 


The data reported here were obtained using six Inlet suppressors used 
singly and togethet on an Avco-Lycomfng YF102 eng1ne» The engine was run 
with a large aft muffler which removed nil aft fan, core, and Jet noise 
and assured that only Inlet noise was being studleu« Each liner section 
had a length to diameter ratio of one-half. The data reported here was 
for the fan blade passage frequency. Additional details concerning the 
tests and the liner construction can be found In reference 7» which also 
contains a large amount of additional data. The parameters for the liners 
for which data Is reported here are found in Table I, The resistance and 
reactance were calculated using? 

n - , l\ 

a”(2 + 1.256 dVd^) 



cot 


2rrfb 
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( 14 ) 


The loci of optimum impedance for the highest speed engine conditions 
for three boundary layer thicknesses, as well as the Impedances of the 
six liners, are shown in Figure 4. These calculations, as well as all 
that follow, are for the blade passage frequency of 4734 Hz. Three bound- 
ary layer thicknesses were considered in Figure 4 to represent an average 
value for each of the liner locations. Liner Bl, being located nearest 
the Inlet Up, is associated with the thinnest boundary layer while liner 
B3, nearest the fan, is located in the thickest boundary layer region. All 
possible propagating modes were considered in generating the optimum Imped- 
ance loci. Mainly axial or well propagating modes with high cutoff ratio 
are located on the left end of the loci while near cutoff (C 1) or mainly 
transverse propagating modes are located at the bottom right. As can be 
seen, liner B3 was designed for near cutoff modes, B2 for Intermediate cut- 
offs ratios, and Bl for well propagating modes. All of the liners were In- 
tended to fall on the loci, bur due to a deficiency In the earlier Impedance 
model used at the time of the design, the liners are all under-damped as 
compared to the target values. This will not affect the comparisons shown 
here. Liner B3 was anticipated to provide the highest attenuation, located 
mostly to the sideline due to the damping and radiation properties of near 
cutoff modes. Liner B2 should provide more modest attenuation somewhat 
nearer to the inlet centerline. Liner Bl should have lower attenuation 
located even more toward the centerline. The above is based upon single 
mode reasoning. With a multitude of modes thefe Is a complicated superposi- 
tion of radiation patterns, and modes at different cutoff ratios will all be 
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attenuated differently. Hence, the real situation will be much more com- 
plicated than described above. However, if the data Is properly taken 
and interpreted, the above qualitative liner behavior can be seen. 

Liners and B5 were included to provide underdamped and overdamped 
versions of liner B3. The variation of far-field directivity attenuation 
with liner resistance can be studied using these three liners. Two sec- 
tions of liner U1 were built to represent somewhat of a baseline, since 
this liner represented our best estimate of a good uniform liner at the 
time the liners were designed. 

Figure 5 shows the experimental results and the theoretical compari- 
sons for a progressive untaping of the Bl, B2, B3 liners. At first, all 
of the liners are taped to establish a hard-wall baseline and to estimate 
the acoustic power versus cutoff ratio distribution as discussed In the 
previous section. When liner B3 Is untaped (located nearest the fan), 
the directivity attenuation shown by the circle symbols is obtained. The 
peak attenuation occurs somewhat less than ip = 70 degrees and is quite 
large (about 9 dB) . The angle shift due to convective effects In the duct 
inlet, as represented by equation (ll), is seen to be essential to obtain 
even qualitative agreement between theory and data. Without this shift, 
the theoretical curve for liner B3 would have peaked at ip <= 90 degrees. 
This would occur since there are many modes near cutoff which are easily 
damped, and according to the zero Mach number flanged duct theory, these 
would radiate near 90 degrees from the Inlet axis. 

With this convective correction as applied in Figure 5 the theoretical 
attenuation directivity agrees quite well with the data. There is an over- 
prediction of about 35 percent at the peaks, which is quite good for this 
initial comparison. The second set of data, labeled B 32 - 83 , is obtained 
by subtracting the attenuation of liner B3 only from that when both B3 and 
B2 are untaped. Thus, liner B2 is working on the modes emerging from 
liner B3, which should be deficient in modes near cutoff. Again, the agree- 
ment with the theory is quite good. The final set of data, labeled B 32 ] - 
B 32 , is again an incremental attenuation due to liner Bl operating on the 
modes remaining after attenuation by liners B3 and B2. Near cutoff and 
moderately cutoff modes should be deficient, leaving only axially or well 
propagating modes. The remaining attenuation is quite small and Is seen 
to be accurately predicted by the theory. The abrupt termination in the 
theoretical curves at about 62 degrees is due to the use of the flanged 
duct radiation model. This model is valid only to 90 degrees for zero 
Mach number and this point is shifted down to near 60 degrees due to the 
convective correction discussed above. Use of an unflanged radiation model 
will alleviate this problem. 
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Although the pr ediicttons are a Ifttle larger than the data In Figure 
5, the qualitative description Is excellent.. It should be noted that no 
corrections have yet been applied for inactive liner area due to fastening 
of the liner perforate to the spool piece flanges. After this has been ac- 
counted for, a simple empirical correction could be applied to provide ex- 
cellent agreement near the peak. This has not yet been done since there 
are first order changes which must yet be made to the approximate theory. 
This will become more apparent In the later discussions. No explanation 
can be made at this time as to why the theory is inconsistent at small 
angles. 

In Figure 5 liner B3 was tested in isolation but liners B2 and B1 were 
tested with a somewhat conditioned power input distribution. That Is, 
other liners had operated upon the power distribution before it entered B2 
or Bl, Data from tests of liners B2 and Bl by themselves are shown In Fig- 
ure 6„ The directivity attenuations are quite different than for the tests 
shown In Figure 5> This is true in particular for liner Bl shown In Figure 
6 (compare to lowest data points in Figure 5) Under the test conditions 
shown in Figure 5, liner Bl had only the well propagating modes to operate 
upon, but when tested in isolation the full range of power for all cutoff 
ratios was available,. The attenuation is much larger in Figure 6 and Is 
biased toward larger angles. This is due to the presence of near cutoff 
modes. 

Notice that the theoretical curves for B2 and B3 both overpredict the 
data in Figure 6, while for Bl the theory under predicts the data. It Is 
this Bl liner result which initiated the near cutoff contour broadening 
study reported In the Appendix. Bef e this modification the theory under- 
predicted the data considerably more than in Figure 6„ The equations (n 
the Appendix represent all of the contour broadening which can be theoret- 
ically Justified. Additional broadening could be empirically inserted, but 
this does not seem to be justified at this point in time. Other theoret- 
ically treatable effects such as modal scattering should be studied first. 

The variation of attenuation directivity with liner resistance is shown 
in Figure 7. The agreement between theory and data is qualitatively correct 
but a greater sensitivity to resistance is exhibited by the theory than by 
the data. The agreement for the two higher resistance liners (B3 and B5) 
could be improved by using an empirical multiplier but then the lower re- 
sistance data would be somewhat underpredicted. 

The attenuation directivity obtained from the progressive untaplng of 
liner Ul is shown in Figure 8. For the first L/D “ l/l| the liner fs seen 
to be quite effective with the theory somewhat underpredicting the peak. 

The theory shows very little change for the second L/0 = 1/A while the 
data shows a considerable drop in liner effectiveness. The final L/D = 1/2 
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of Hnef Vopfng shows a comparison botwean theory and data more like the 
higher attenuation curves of previous figures, with the theory somewhat 
overpradicttng the data. 

An attempt was also made to predict the attenuation directivity of 
multiple pure tones at takeoff engine speed. Problems are Immediately 
apparent, as Illustrated In the narrow-band far-fleld directivity In Fig- 
ure 9 for the twelfth shaft harmonic. The only simple assumption to make 
as to the ntodal content of this noise Is that It exists as a 12-lobe pat- 
tern (m •• 12), The theoretical directivity of the 12-lobe mode, first 
radial, Is shown on the figure and the agreement looks plausible for high 
angles. However, the noise evident between 10 and ItO degrees Indicates 
the presence of one or nrare other modes which are not so easily deduced* 

The sound attenuation obtained experimentally with a liner designed to 
remove multiple pure tones (ref, 7) reached about 18 decibels for a liner 
length of L/D « l/l». However, the theoretical attenuation using equal 
acoustic Dower per mode (which Is not considered a good assumption for 
MPTs) reached only about 5 decibels, With the assumption that only the 
m *• 12 mode Is present, the theoretical attenuation Is also only about 
5 decibels. It Is obvious that In Its present form, this method cannot 
be used for multiple pure tones where only one or a few modes are dominant. 
It Is clear that modal scattering plays an Important role here, and the 
following example emphasizes this point. The 12-lobe pattern has ^ » 

1.0l<2 In the hard-wall duct. At the liner optimum Impedance for this mode, 
the lowest radial order has ^ « 0.807 with an attenuation of 2.2 times 
that predicted by equation (2). The higher radial modes, while not Investi- 
gated, may have even higher attenuation. Thus, It Is seen that modal scat- 
tering must be Incorporated Into this program If It Is to be used for multi- 
ple pure tones. Of course, exact propagation analysis can also be used In 
this case Involving only a few Incident modes. 


CONCLUDING REMARKS 


This paper has presented the Initial comparison between experimental 
suppression directivity data and the theoretical prediction?* using the 
simplified liner evaluation program based upon mode cutoff ratio. It 
should again be pointed out that the prediction of attenuation directivity 
Is an extremely demanding objective. This prediction Involves noise source 
considerations, duct propagation, and far-fleld radiation from the duct 
termination. In spite of this complexity, the theoretical model provided 
reasonable quantitative and excellent qualitative agreement with the exper- 
imental data. The results of the theory presented here did not have any 
empirical corrections applied. Empirical correction Is a path that will be 
purs led. There are many elements In this program which may be Individually 
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corrected and tested as experimental evidence points to the need. The 
theoretical model can also be improved. Inclusion of modal scattering 
at the liner input may be of benefit. This has been shown to be essen- 
tial for possible single mode propagation such as multiple pure tones. 

The addition of an unflanged duct radiation model will also provide an 
improvement. Also, the potential change? <n ;he radiation pattern due 
to modal scattering through duct area ch.mgei such as ti e throat and the 
Inlet lip must be Investigated. Several if the above suggested additions 
to the theory, us well as a more careful evaluation of the sou'^'ce power 
distribution, may provide a better agreement between theoretical and ex- 
perimental far-field directivity attenuation. 
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APPENDIX “ OFF OPTIMUM ATTENUATION 


The method for calculating liner attenuation used In this report does 
not Involve exact modal propagation calculations, but Instead Is based 
upon correlations derived from these calculations. The correlations In- 
volve the values of the maximum possible attenuation ( AdBn,/L/D) as expres- 
sed In equation ( 2 ) In the text and of the optimum Impedance components 
(0m» ^m) based upon the equations in reference If. Since In general the 
Impedance of Interest In a calculation Is not the optimum Impedance, a 
method must be provided to calculate attenuation at any arbitrary impedance. 
Such a method was previously reported In reference 10, which provided con- 
stant damping circular contours In the Impedance plane to approximate the 
exact contours. An example of these circular contours was shown In Fig- 
ure 1. This correlation was based upon well propagating modes (large C) 
and did not consider the effects of a boundary layer near the wall, 

Comparisons of theoretical and experimental far-fleld attenuation direc- 
tivities as presented In this paper Indicated that the correlation did not 
provide sufficiently broad attenuation contours for modes near cutoff. A 
closer examination of the exact attenuation contours for modes near cutoff 
showed that thisvvas indeed the case. Also discovered in this reexamination 
of exact theory is that for well propagating modes, when the boundary layer 
caused a large reduction of optimum resistance (see ref. k) , the entire sys- 
tem of attenuation contours was compressed. 

The new correlation, an example of which is shown in Figure 2, accounts 
for near cutoff broadening and compression due to the boundary lay-'i*. The 
basic equation for the circular contours Is, 

(0 - 0 j-)^ + (x - Xc)^ (Al) 


where the coordinates for the cent.er of the circle are, 


0 . 


0 + FY^(Y-l)(gV 2 . 1 ) / ®m' 


m 


(Y^+7) 


W~ 

, mo/ 


(A2) 


and 


V == Y 
'‘c "^m 


(A3) 
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The radius of the circular contour Is, 




jhIHIL X FY^h-IHeVz.l) 

(Art 7 ^ 



M) 


where 


F 


8.7 


AdBm 

TTd 


(l+Mp)^ 


and 


Y 


AdBm 
“ AdB 


(A5) 


(A6) 


which Is the ratio of maximum possible attenuation to the actual attenua- 
tion on the contour. 

The quantity 0m/0mo represents the effect of contour compression due 
to boundary layer refraction. The optimum resistance with a boundary layer 
(0^,) , without a boundary layer (0mo) optimum reactance (Xj,)) are all 

found using the correlations of reference 4. Equation (2) In the text pro- 
vides the value of AdB^^L/D. 

Equations (A2) , (A3), and (A4) can be used to construct the circular 
contour system as Illustrated In Figure 2. In the actual attenuation pro- 
gram used In this paper, equation (Al), using equations (A3-) , (A4) , and 
(A5) , was solved by Iteration to find the attenuation at any arbitrary Im- 
pedance. All quantities In the equations are known except AdB In equa- 
tion (A6) for Y* To start *,he Iteration y Is set equal to one and then 
Increased until equation (Al) Is satisfied. With y thus determined, the 
desired attenuation (AdB) at the arbitrary Impedance Is found from equation 
(A6). 
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Figure 2. - Comparison of exact and approximate equal damping con- 
tours lor a near cutoff mode, m • 29, |i ■ 3, ^ ■ 1,0626, • 14. 2 
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Figure 3. • Optimum Impedance Iqcus with equal mode number Intervals. 
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Figure 4. - Calculated optimum Impedance loci and wall Impedances for three suppressors 
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